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Measurement of collisions between rubidium atoms and optically dark rubidium ions 

in trapped mixtures 
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We measure the collision rate coefficient between laser cooled Rubidium (Rb) atoms in a magneto- 
optical trap (MOT) and optically dark Rb + ions in an overlapping Paul trap. In such a mixture, 
the ions are created from the MOT atoms and allowed to accumulate in the ion trap, which results 
in a significant reduction in the number of steady state MOT atoms. A theoretical rate equation 
model is developed to describe the evolution of the MOT atom number, due to ionization and 
ion-atom collision, and derive an expression for the ion-atom collision rate coefficient. The loss of 
MOT atoms is studied systematically, by sequentially switching on the various mechanisms in the 
experiment. Combining the measurements with the model allows the direct determination of the 
ion-atom collision rate coefficient. Finally the scope of the experimental technique developed here 
is discussed. 

PACS numbers: 34.50.-s,52.20.Hv, 37.90. +j,37.10.Ty 



I. INTRODUCTION 

Trapped ion and atom mixtures enable us to inves- 
tigate exciting physics ranging from, a single trapped 
ion inside a Bose-Einstcin Condensate (BEC) [TJ [5] and 
Magneto-Optical Trap (MOT) atoms superimposed with 
laser cooled ions [3H5] , to mixtures where the ions are 
only cooled in collision with the atoms [SHE]- Elastic 
and charge exchange processes in these mixed, trapped 
system have been demonstrated and measured [H [5] . 
Even before the existence of such mixed traps, a variety 
of rate coefficients and cross-sections for several ion-atom 
processes have been measured [TU]. The most prevalent 
method for the determination of the ion-neutral collision 
rate using trapped ions, measures the loss of ions from the 
trap [TUHH]- For optically dark ions, this loss is measured 
by sweeping the ion resonance frequency or by extracting 
the trapped ions to charge particle detectors [ID]. How- 
ever, in mixed ion-atom trapping experiments in which 
the ions are optically dark, the cold trapped atoms in- 
teracting with ions can be directly measured in situ, by 
atomic fluorescence, without affecting the trapped ions, 
fn this article, we develop a method to measure ion-atom 
collision rate coefficient in such a hybrid system using 
atom fluorescence as a probe. 

fn the experiment reported here, we trap the opti- 
cally dark Rubidium ( 85 Rb + ) ions, derived from the laser 
cooled 85 Rb atoms in a MOT. The ions are created and 
accumulated at the center of the ion trap, which overlaps 
with the MOT. Such a system evolves into a mixture with 
constant numbers of trapped atoms and ions in steady 
state. The ionization of the MOT atoms and the sub- 
sequent trapped ion interaction with the atoms results 
in a depletion of the atoms from the MOT. The conse- 
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quent drop in MOT fluorescence is utilized to develop 
the framework for the measurement of the collision rate 
coefficient for ion-atom interactions, as described below. 
In what follows, we construct the method for the 
determination of the collision rate coefficient. We first 
present a brief description of the experimental setup, fol- 
lowed by a simple, analytical rate equation model, where 
we systematically consider loading and loss processes of 
the MOT atoms. Here we model in a calibrated manner 
the change in the trapped atom number when each ion 
related process is turned on. This culminates with the 
simultaneous trapping of ions and atoms in steady state. 
The experimental data is then presented in a sequence 
consistent with the processes modeled. Combining the 
measurements with the derived analytical expressions al- 
lows the determination of the rate of atom loss due to 
ion-atom interactions, and hence the collision rate coef- 
ficient for ions and atoms. We end with a discussion of 
some features of such measurements. 



II. EXPERIMENTAL SETUP 

The experimental arrangement (Fig. nl), techniques 
and characterization have been described in detail in pre- 
vious work |12j . Here we include a brief description for 
making the experimental discussion comprehensible. Rb 
atom vapor is created from a dispenser and 85 Rb atoms 
are laser cooled and trapped in a six-beam MOT of size = 
0.35 mm (standard deviation) and temperature = 166(± 
28) fiK. The cooling laser is red detuned by S c — 2tt x 10 
MHz with respect to the 5s 1/2 (F = 3) -B- 5p 3/2 (F' = 4) 
transition, and the laser power is 3.5 mW, per beam. A 
magnetic gradient field of 12 Gauss/cm is utilized for the 
MOT. The intersection point of the MOT laser intensity 
maxima and magnetic field zero coincides with the center 
of a linear Paul trap, built around the capture volume of 
the MOT. This ensures that the trapped ions and laser 
cooled atoms are optimally overlapped [SJ [T^] . Ions are 
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FIG. 1: (color online) The schematic diagram of the exper- 
imental setup is illustrated. The MOT and the linear Paul 
trap are overlapped as shown, with the cooling lasers illus- 
trated in red, and the MOT located at the intersection point 
represented by a red sphere. The ionizing blue light is illus- 
trated, without its source. The channel electron multiplier 
(CEM) is used to detect the trapped ions, by appropriately 
switching the voltage on the hollow end cap electrode closer to 
the CEM. The direction for the MOT fluorescence detection 
is indicated, without illustrating the spatial filtering arrange- 
ment. 



III. MODEL FOR THE LOSS OF TRAPPED 
ATOMS 

A. The MOT rate equation 

We now present the rate equation framework which 
describes the evolution of the ion atom system as we add 
the various ingredients of the above experiment one at a 
time. The total number of atoms in the MOT (Nmot) 
at any given time, is Nmot = N g +N e , where N g and N e 
are the atom numbers in the ground and excited states 
respectively. The time dependence of Nmot can be writ- 
ten as 



dN M QT 

dt 



= L- JmlN J 



MOT, 



(1) 



with L, the loading rate of the atoms from the back- 
ground vapor into the MOT and j m i, the loss rate of 
the MOT atoms, which is primarily collisional [IS]. The 
time dependent solution of eqn (fTl) with initial condition 
N M ot(0) = is 



N MOT (t) = N (l - e~™*), N Q = 



L 

1ml 



(2) 



where No is the steady state atom number in the MOT, 
in the limit of large time. 



B. Fraction of excited atoms 



created from cold atoms using a blue light source (BLS) 
which along with the cooling light results in two photon 
ionization of the Rb MOT atoms. The ion trap is made 
up of four parallel rods in the quadrupolc configuration, 
where a time varying radio-frequency (RF) voltage is ap- 
plied to one diagonal set, and an opposite phase RF volt- 
age to the complementary diagonal set. This generates a 
two dimensional quadrupole field for the radial trapping 
of ions. Two end cap ring electrodes at a positive dc volt- 
age confine the positive ions along the trap axis. Atoms 
arc detected by observing the fluorescence from the MOT 
on a photo-multiplier tube (PMT) and imaged using two 
CCD cameras. The MOT fluorescence is measured by the 
PMT, where a spatial filtering arrangement is used to re- 
duce background light signal. The ions are optically dark 
and are therefore detected either indirectly by a change 
in atomic fluorescence of the MOT or destructively by 
extraction onto a channel electron multiplier (CEM) and 
measuring the ion counts, either by pulse counting or 
by measurement of the proportional extracted ion signal. 
The ion extraction to the CEM is affected by switching 
the ion trap end cap electrode from +80 V to —5 V, while 
the RF field is still on. 



The fraction of atoms in the excited state f e is de- 
termined from the photon scattering rate from an atom 
and can be written as [T3] 



fe = 



N, 



(n/rf 



N g + N e i + 4(<s c /r) 2 + 2(ft/r) 2; 



(3) 



where, fl = 2-zr x 16.7 MHz is the Rabi frequency and 
r = 2-/T x 6 MHz is the decay rate of the excited atomic 
state. 



C. Photo ionization of atoms 

The ionization of atoms in our experiment is accom- 
plished by excitation of the atoms in the 5p 3 / 2 state to the 
ionization continuum using the BLS. The central wave- 
length the BLS is at X P i = 456 nm, which is sufficient 
to ionize the excited Rb atoms 16 . The intensity of 
the BLS light incident on the MOT is denoted by J p j. 
For photo ionization energy E p i — hc/X P i, the ioniza- 
tion cross-section is denoted by a p i. The photon energy 
which is in excess of the ionization energy is carried away 
by the ejected electron, with the resulting ion suffering 
the recoil. 

The loading rate of the MOT remains unchanged 
when the BLS is switched on for photo ionization, as the 



population of the background atoms remains unchanged 
with operation of the BLS. The maximum intensity of 
the BLS used is 1.33 mW/cm 2 . The rate equation for 
the MOT, when the BLS is ON is modified to 



(INmot 
dt 



L - ImlNMOT — IpiNMOT, 



(4) 



where a new loss rate -f P i adds to the losses inherent to 
the operation of the isolated MOT to give a total loss rate 
of -f t = j m i + -f pi [T71 [T5] . Here we assume that since the 
ion trap is switched off during photo ionization, the ions 
and electrons created from the MOT atoms immediately 
leave the system and play no further role in determining 
the number of trapped atoms. The ionization rate of 
the excited Rb atoms by radiation from the BLS can be 
written as IT71H81 



l P i 



E - 



-*-pi S-'-'j 



pii 



(5) 



where C, — ( a "\^. ) • The corresponding time depen- 
dence of Nmot with initial condition Nmot(0) = and 
the inclusion of the ionizing process can then be written 



N MOT (t) = —{l-e 

It 



-itt\ 



(6) 



D. Ion- Atom interaction loss 



Once the ion trap is switched ON, then a fraction 
of the created ions load into the ion trap which overlaps 
the MOT. Since the ion and atom traps are optimally 
overlapped, the ions created from the MOT have negli- 
gible velocity upon creation. These ions gain in kinetic 
energy as they evolve in the field of the ion trap. As the 
number of ions loaded in the ion trap grows, the spa- 
tial distribution of the trapped ions increases outward 
from the trap centers. Since the MOT volume is much 
smaller than the volume of the trapped ion spatial distri- 
bution, therefore long before the total number of trapped 
ions equilibrates, the number of ions at any instant over- 
lapped with the MOT stabilizes. Therefore the collisions 
rate between the trapped ions and the MOT atoms also 
stabilizes rapidly, allowing the definition of a time in- 
dependent binary ion-atom interaction rate, -fi a for the 
constituents of the ion-atom mixture. Since atoms are 
weakly trapped with respect to ions, they can gain suf- 
ficient energy in collision with trapped ions to exit the 
MOT. This is a new loss channel for the atoms trapped 
in a MOT, which can be written as, 



dN 



MOT 

~dt 



{lm.1 + Ipi + lia)N M OT- 



(7) 



The resulting time evolution for Nmot with initial con- 
dition Nmot(0) = therefore becomes, 



N 



MOT 



(t) 



Itot 



(1 



-ftott 



(8) 



where jtat — It+lia is the total loss rate of MOT atoms, 
when the ion and atom traps are operated simultaneously 
as described in the experimental section above. We note 
that while the ion density overlapping with MOT volume 
stabilizes rapidly, the number and the velocity distribu- 
tion of the trapped ions evolve continuously till the ion 
trap is filled to its limit. Thus the rate of ion loss from 
the ion trap and the evolution of Nmot are linked. 



E. Number of trapped ions 

The ion trap is loaded to the minimum of the secu- 
lar trap potential, by photoionization of the overlapping 
MOT atoms. The number of trapped ions as a function 
of time, for a given photo- ionizing intensity I p i, can be 
described as 



Ni(t) = Nmot C — (1 

7 



--ft 



), 



(9) 



where 7 is the ion trap loss rate. Since ions have much 
higher velocities than the MOT capture velocity for the 
atoms, ion-atom collisions will lead to loss of atoms from 
the MOT. This allows the number of MOT atoms lost 
due to trapped ions N M q T , to be written as 



N M s 3 T (t)=N(l-e-^% 



(10) 



where N is saturation number of atom loss due to ion- 
atom interaction, and -fi a is the loss rate of the MOT 
atoms due to ion-atom interactions. The above form of 



the atom loss will be validated later (eqns. (21 1,(22) and 
Fig. pi), when the experiment and results are presented. 
The atom loss rate from the MOT, due to ion atom in- 
teractions, is proportional to the number of trapp ed ions 



(N M s £ T (t) ex Ni(t)), and so from eqns. (|9| and (|10|) we 



conclude that 7 = -fi a - Therefore, number of trapped 
ions as a function of t is expressed as, 



Nj{t) = Nmot ( ^ (1 

lia 



-lia,t 



(11) 



The number of trapped ions also depends on I p i but can- 
not increase indefinitely with I p i due to finite trap depth 
and volume. Therefore an intensity loss coefficient of 
trapped ions must be introduced. In the same spirit as 
the above discussion, to construct a general form of A/, 
we write 



dNjitJpi) 



dl, 



N 



1 



MOT 



C—0.-e-*"* t )-KN I (t,I pi ), 



pi 



(12) 



where k is the intensity loss coefficient and has the unit 
of inverse intensity. The solution of the above equation 
is, then 



1 



N I (t,I pi )=N MO Tt (1 



-liat 



JiaK 



)(1 



-K,l pi \ 



(13) 



The equation can be fully converted into a time depen- 
dent function of trapped ions when I p i — > oo and also 
fully converted into an I p i dependent function of trapped 
ions when t — > oo. The number of ions trapped as a func- 



tion of Ipi for t — > oo is written as 



N r (I pi ) = JV?(1 - e~ K H 



(14) 



The maximum trapped ion number, Nj when t — > oo and 
I p i — >• oo is expressed as, 



N° 



Nmot C 



liaK 



(15) 



Since eqns. (14 1 and (151 apply to the ions system in 



the limit t — > oo, it is therefore reasonable to expect an 
average motional energy (Ei), for each ion. 



F. Determination of ion-atom collision rate 

The ions trapped as described above are assumed 
to have a speed distribution f{v), consistent with (Ej). 
The laser cooled atoms on the other hand are expected 
to exhibit a temperature Ta ~ 100 pK. In this situation 
{Ei) >■ ksTA and therefore all of the velocity of the ion- 
atom collision can be assumed to be possessed by the 
ion, in the lab frame of reference (LFoR). Similarly, as is 
shown later, Nmot 3> Nj as the ion density pi, is far 
less that the atom density pa in such mixtures despite 
the ion trap volume, Vit ^> Vmot, the MOT volume in 
our experiment. 

The total ion-atom collision cross-section, tr tot is en- 
ergy (E oc v 2 ) dependent and is the sum of the elastic, 
a e i, and resonant charge exchange, a cx . a e i oc 1/E 1 ' 3 
over all the energies and a cx oc 1/E 1 / 2 in the low energy 
regime and a cx oc (a In E — b) 2 at high collision energies 
(a and b are constants) [19]. Given f(v) for the ions, 
the determination of Otot rests on assumptions made for 
f(v). In this case the experimentally meaningful quantity 
to consider is the ion- atom rate coefficient, 



h- — 



<?tot V 



f{v)dv 



{crtot v ; 



(16) 



which represents the velocity averaged product of <7 tot t>. 
A single ion with velocity Vq and a corresponding ion- 
atom cross-section of ot will collide with MOT atoms of 
density pa at a rate given by 



z = Vtot v o PA- 



(17) 



For the Ni-m ions that overlap the MOT at any given 
time, with a speed distribution f{v), the total ion-atom 
collision rate then becomes, 

zNi-m = {otot v) paNi-m = k ia Pa Ni- M - (18) 

For a given value of I p i, the number of ions in the ion 
trap volume Vit, overlapping the MOT volume Vmot, 
is given by, 



N 



I-M 



jy?(i 



^ — kI~ 



') 



Vmot 
V, 



it 



(19) 



Since ion velocities are large, and the MOT capture 
velocity for the atoms much smaller, ion-atom collisions 
cause the MOT atoms to eject, allowing us to equate the 
total ion-atom collision rate in eqn. ( fTS] ) to 7i a ^A/OT, the 
atom loss rate due to ion-atom collisions. After substi- 



tuting eqn. (191 into eqn. (18), the resulting expression 
for "fin then becomes, 



7i 



= N°k la 



V, 



IT 



^ — kI 



'"), 



(20) 



which defines the relation between the experimentally 
measurable ji a and the rate coefficient for ion-atom colli- 
sions. In the experimental results that follow, we demon- 
strate that the above rate equation formulation describes 
adequately the dependence of the various loss rates with 
Ipi, the BLS intensity. 



IV. RESULTS 

A. Experimental Sequence 

We now describe the experiment which allows us to 
validate the above rate equation analysis for ion-atom 
interaction. The basic time sequence instrumental for 
results below is shown in Fig. [2] Here a MOT containing 
f=a 1.80(±0.06) x 10 6 atoms is loaded to saturation in 40 
s (Ti to T2), and the change in the number of atoms 
in the MOT is measured under two conditions from T2 
to T3 at different intensities of the ionizing light (BLS). 
When the BLS ionizes the cooled 85 Rb atoms from the 
MOT and these ions and electrons leave the MOT volume 
immediately as there is no confinement for them. This 
results in a depletion of the steady state atom number in 
the MOT and therefore allows the determination of the 
ionization rate for the excited atoms. Using eqn. pi, data 
from Ti to T2 are fitted to determined 7 m ;. Solution for 
eqn. Q with appropriate initial conditions, N(t = 0) = 
N and N(t -> 00) = N a - N t is 



N 



MOT 



{t)=N -N t (l 



-~ttt 



(21) 



where Nt is the atom loss due to 74 during photo ion- 
ization. 7 t is determined by fitting data from T 2 to T 3 
in Fig. [2ja) with eqn. Q. In Fig. [2|b) we switch ON 








10 



20 30 40 

time (s) 

T 2 



50 



60 




10 



20 



30 40 
time (s) 



50 



60 



FIG. 2: (color online) The temporal sequence for the exper- 
iment is illustrated. In both panels, first the MOT is loaded 
for 40 s (from Ti and T2) to saturation, where it contains 
~ 1.8 x 10 6 atoms. Beyond T2, the photo ionization light 
is switched on. Panels (a) and (b) illustrate the evolution of 
the trapped atoms when the ion trap is OFF and with ion 
trap ON, respectively. The number of trapped atoms reduces 
in both cases. However, when the ions and the atoms are 
simultaneously trapped, the atom loss is much more signifi- 
cant. The plots in (a) and (b) correspond to BLS intensity 
I P i — 1.33 mW/cm 2 . The various loss rates discussed in the 
text are determined from such data as a function of photo- 
ionizing intensity. N t and Ntot are the atom losses due to 
7t and ytot- The lines show the fits to the data according to 
eqns. (211 and ( 22 I repectively. 



both the ionizing light (BLS) and the ion trap simulta- 
neously. In this case the ionization process is unaffected, 
but some fraction of the ions created reside within the 
ion trap, which is well overlapped with the atom trap. 
Here we observe that the number of trapped atoms re- 
duces significantly over the ionization loss. This loss is 
attributed to ion-atom interactions. Solution for eqn. flh 
is obtained in a similar way as photo ionization case and 
given below. 



N. 



MOT 



= N - N tot {l - e-**-'*) 



(22) 



where N tot is the atom loss due to j tot when both atoms 
and ions are simultaneously trapped. 




BLS intensity (mW/cm ) 



FIG. 3: (color online) The loss rate of atoms from the MOT 
as a function of photo ionization light intensity, I pi , with- 
out the ion trap in Fig. El (a). As I p i increases, the loss of 
atoms from the MOT increases as does the rate of loss ■jpi. 
Since this atoms loss is happening on top of the normal loss 
from the MOT, the experimentally measured quantity here 
is 7 t . As can be seen from the text, a linear relation is ex- 
pected between -yt and I p i, which is measured with a slope of 
0.00225(±0.0002) m 2 /J and the 74 axis intercept of « 0.14 s" 1 
which is the MOT loss rate without the ionizing light (BLS). 
Error on the slope from the fit is calculated from confidence 
level analysis. 



B. Ionization without trapping 

When the BLS is switched ON at T2, and the ion 
trap is not operational, the loss rate of atoms from the 
MOT increases due to ionization of the atoms, which 
is shown in Fig. [3] The measured loss rate coefficient 
It = "Imi + Ipi can be written using eqn. (15J) as, 



It = 1ml + (I P 



(23) 



From Fig. [3] it is clear that 74 has a linear dependence on 
the intensity of the BLS, i™. This is because although the 
ionization process is a two photon process, the operation 
of the MOT ensures that a constant fraction of atoms 
are present in the excited state, making the ionization 
an effective single photon process. The slope of 7 t , £ is 
determined to be 0.00225(±0.0002) m 2 /J from the data. 
The value of j t for I pi = is j ml = 0.14(±0.011) s _1 . 



C. Ionization and ion trapping 

For the case when ions are held in the ion trap after 
T2 in Fig. pjb), loss rate coefficient ^ tot varies with a 
nonlinear behavior as shown in Fig. |4j The change in -f t ot 
is fit by a combination of the linear expression in eqn. ( 23 ) 
and the exponential form of eqn 



( |20[ ). To study the loss 
rate due to ion-atom interaction, ji a is plotted separately, 
as ji a — 7tot~7t) m Fig- 00 The exponential nature of the 
7i a as a function of I p i is evident thus validating the rate 
equation for the ion-atom interaction. From the fit to the 




BLS intensity (mW/cm ) 

FIG. 4: (color online) The loss rate of atoms from the MOT 
as a function of ionization light intensity, I p i, with the ion 
trap in Fig. p] (b). As I P i increases, the loss of atoms occurs 
on top of the normal loss from the MOT and the ionization 
loss, and so the experimentally measured quantity here is "ftot- 
The additional loss of atoms over that in Fig. 13] is attributed 
to the ion-atom interactions for which the loss rate is ^i a . "/tot 
has exponential and linear components in I p i as is observed. 
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FIG. 5: (color online) A subtraction of the data in Fig.^lfrom 
that in Fig. [4] is plotted to isolate the ion-atom interaction 
term. This loss rate of the MOT atoms due to ion atom 
interaction, 7; a is plotted against BLS intensity, I P i. 7i 2 is 
fitted with a single exponential function of the form, eqn. ( 20 1 
and the coefficients of the fit are utilized to characterize the 
interaction term, as discussed in the text. 



data eqn. (20) becomes, 7i a = 0.218(1 — exp[— 0.126 I p 



Identifying k — 0.126(±0.031) in inverse intensity units, 
i.e. m 2 /W and 



Vit 



= 0.218(±0.029) 



(24) 



with units s _1 , we can constrain the product of the 
trapped ion density and the ion-atom collision rate coef- 
ficient. 



D. Direct ion detection 

Alongside the measurements on the MOT discussed 
above, the trapped ions though optically dark in the 
present experiment, can be measured using a CEM as de- 
scribed in the experimental arrangement. At the end of 
each ion- atom experiment, the trapped ions are extracted 
into the CEM, by switching the end cap voltage appro- 
priately. However, since this experiment fills the ion trap 
to its capacity, severe pile-up results due to overlapping 
arrival times of the ions onto the detector. For a CEM 
cone voltage of -2100 V, the extracted ions from the trap 
as a function of I p i is illustrated in Fig. [6J A single ex- 
ponential dependence of the ion numbers detected with 
BLS intensity can be expected and a dependence of that 
nature is seen in the data. However there is detector sat- 
uration to contend with and the present measurement of 
ion numbers is limited by the detection. The direct mea- 
surement of the CEM signal is therefore presented for the 
sake of experimental completeness rather than with the 
motive of utilizing it to determine the ion-atom interac- 
tion rates. For large ion numbers this is likely to remain 
a significant problem with present day ion detectors. Al- 



ternative schemes for ion detections by electrode pickup, 
Faraday cup detection, etc. [11], may be useful for this 
problem and should be explored. 



E. Collision rate coefficient determination 



We now demonstrate the use of eqn. (24), to arrive 



at the ion-atom collision rate in the present experiment. 
The values of all the quantities used in the determina- 
tion of the collision rate, with the standard deviation 
error in the accompanying parenthesis are provided in 
the text below. The product of the ion density and ki a 



is determined from the fit to eqn. Q as 0.218(±0.029) 
/s, so only the determination of the number of trapped 
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FIG. 6: (color online) The change in the detected ion counts 
vs. BLS intensity is shown. Since the number of trapped ions 
is very large, severe pileup of the ion signal results and the 
ion signal is seen to saturate faster than measured from the 
atom signal. The fit is generated using eqn. fl4|). 



TABLE I: The key values and errors for the quantities which are used in the determination of ki, 

Method 

fit from the experimental data, in Fig. |3| 
fluorescence measurement 

derived from secular frequencies and trap depth 
refer eq. 



Parameter 


Value (± error) 


Units 


c 


0.00225(±0.0002) 


m 2 /J 


Nmot 


1.80(±0.06) x 10 6 


number 


Vit 


8.33(±0.83) x fO" 8 


m 


N° 


1.48(±0.43) x fO 5 


number 


lia{Ipi -> OO) 


0.218(±0.029) 


s^ 1 


K 


0.126(±0.031) 


m 2 /W 




1.23(±0.42) x ffr 13 


m 3 /s 



15 



fit from the experimental data, in Fig. 5 
fit from the experimental data, in Fig. 5" 
calculated from the various quantities given above 



ions and the volume of trapping for the ions is required 
for computation of the collision rate coefficient. Keeping 
the compromised CEM data in mind, k — 0.126(±0.031) 
m 2 /W is determined from the experimental data for atom 
loss Fig. [5j rather than from the data in Fig. [6] Using 
eqn. (151, the maximum number of ions which can be 



accumulated in the ion trap volume is calculated to be 
N? = 1.48(±0.43) x 10 5 . 

The volume within which the ions are trapped, Vit, 
is determined for the present trap parameters by a com- 
bination of trap loss and trap secular frequency measure- 
ments. A Monte-Carlo analysis of the time for the loss 
of ions from the ion trap [B] gives the secular trap depth 
for trapped ions to be « 0.8 eV. The secular frequen- 
cies for the ion trap are u> x — u y — 2ix x 135 kHz, and 
bj z = 2tt x 27 kHz, which allow us to calculate the trap ex- 
tent in each direction from the relation miuj^r 2 /2 = 0.8 
eV, where r € xq,i/o,Zq is the extremal displacement in 
each direction and mi is the mass of the ion. The trap- 
ping volume defined by these dimensions for a single ion 
then is Vit = 8.3(±0.83) x 10" 8 m 3 . 

The collision rate is determined by substituting the 
mean values of the quantities above in eqn. (24), to be 



(1.23 ±0.42) x 10~ 13 m 6 /s. 



(25) 



The standard deviation error is <Jk la = vlZi °f = 
0.42 x 10 -13 m 3 /s, where cr^'s are the individual errors 
for each contributing parameter, shown in Table. IT] This 
demonstrates the technique's ability to arrive at the rate 
coefficients for ion atom processes, even when the ions 
are not directly detected. The above rate coefficient is 
determined when the fraction of atoms in the excited 
state for f e w 0.28, from eqn. [3] The measured rate 
coefficient incorporates the elastic and resonant charge 
exchange collision from the ground and excited state of 
the atoms. 



Theoretical estimate of the collision rate 
coefficient 



To estimate the collision rate coefficient in eqn. ( 16 ) 



we adopt the analytical form of total cross section as 
a function of collision velocity from [THl HO] , using the 
ground state C4 value for Rb atoms, 5.26xl0 _39 Cm 2 /V 



[13] . The speed distribution of trapped ions is deter- 
mined using Monte-Carlo analysis and molecular dynam- 
ics simulations [Hj . Molecular dynamics simulation shows 
that the ion trap can stably trap ions with Maxwell 
Boltzman(MB) distribution of temperature f=s 1000K. In 
Monte Carlo analysis, a non interacting distribution of 
ions are evolved in the absence of cold atoms in the ion 
trap potential to estimate the maximum trappable secu- 
lar energy, which is found to be « 0.8 eV, corresponding 
to speed of 1360m/s. This value matches the tail of the 
MB distribution of ions in the trap. Since Vit >• Vmot, 
and the two traps are well centered, in the overlap re- 
gion, the micro motion velocity is small. Therefore, to 
estimate the speed distribution of the ions, only consid- 
ering the secular velocity is a good approximation. By 
performing the integration in eqn. ( 16 ) with these quanti- 
ties, we calculate the rate coefficient as 9.4xl0 -14 m 3 /s, 
when all the Rb atoms are in the ground state. Because 
atoms in the MOT are constantly pumped to the excited 
state by the cooling laser, ions invariably collide with 
a fraction of the MOT atoms, which are in the excited 
state atoms. In the experiment, we have 28% of atom 
population in the excited state (/ e = 0.28). The scalar 
polarizability for a Rb atom in the excited 5p 3 / 2 state 
is 14.15xlO~ 39 Cm 2 /V 13 , which is larger than that in 
the ground state. Incorporating this fraction of atoms in 
the excited state in the calculation, we compute the rate 
coefficient as ki a = 1.12 x 10~ 13 m 3 /s, which agrees well 
with the experimental value. 



V. DISCUSSION 

The above technique allows the detection of colli- 
sional processes between trapped ions and atoms. In the 
present case, this is successful even for optically dark 
ions. At its core it relies on the measured changes in 
the fluorescence of the MOT atoms in the presence of 
the ions. The present work relies on two important and 
reasonable premises, which need to be emphasized. 

The first is that ion-atom collisions are two body 
processes. This is reasonable because of several condi- 
tions that exist in the experiment. The atoms in the 
MOT are non-interacting to a good approximation. Be- 
cause the ions are hot, they have large velocities and 



in this regime the binary ion-atom cross section is small. 
This is enough to make the presence of another ion in the 
vicinity of the collision, such that it affects the details of 
the ion-atom collision, highly improbable. It should be 
noted though that the experiment and the model are in 
good agreement, and since the model only incorporates 
two particle ion-atom collision, it is fair to argue that 
any more complex process is absorbed within the errors 
quoted in the present result. 

The second premise is that ion-atom collisions knock 
out a MOT atom. In the experiment, the only cool- 
ing mechanism for the ions is by collision with the cold 
atoms [5] . The ions are also subject to significant contin- 
uous RF heating as their numbers are large and Vit ^> 
Vmot, while the cooling is most effective only at the cen- 
ter of the MOT. Thus the ions are quite hot and there- 
fore possess large velocities. Under such circumstances, 
even a glancing ion-atom collision, whether elastic or res- 
onant charge exchange, will transfer sufficient energy to 
the laser cooled atom such that it exceeds the capture ve- 
locity of the MOT and gets ejected from the MOT. This 
then forms the basis for equating the MOT loss rate with 
the ion-atom collision rate as described above. 

Since the majority of ion-atom combination exper- 
iments use visible, laser cooled ions, we provide a brief 
discussion of some major differences. In such cases, since 
ion temperatures are easily in the mK levels, the ion- 
atom cross sections are much higher. The ions then 
would be crystallized and have a velocity distribution 
that would be well characterized. How atoms interact 
with such crystallized ions and what part the long range 
order of the ions plays in the measurement of the binary 
ion-atom interaction needs to be carefully understood. A 
possible problem is that an ion-atom collision at these en- 



ergies may not result in the ejection of the colliding atom 
from the MOT and would therefore require eq. ( 18 ) to be 



written with a proportionality constant. However if such 
an equation can be written, the cross-section could be 
directly determined because of the well characterized ve- 
locity and density distributions of the ions. Obviously, 
the rate equations constructed above would need to be 
modified according to the specifics of such a system. 



VI. CONCLUSION 

In a trapped ion and atom mixture, where the atoms 
are laser cooled and contained in a MOT and the ions are 
trapped within a Paul trap, we have developed a tech- 
nique which measures the rate coefficient of ion atom 
collisions, even when the ions are optically dark. The 
rate equation formalism has been systematically devel- 
oped for this experimental system. The experimental re- 
sults are then fit to the rate equation model, which is seen 
to provide a consistent and adequate description for the 
measurements. Analysis of the experimentally measured 
quantities with the model allows the determination of the 
trapped Rb + ions and Rb atom collision rate coefficient. 
The value of the rate coefficient compares well with the 
theoretical estimation and validates the technique. 
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